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ABSTRACT 

Context. The clustering properties of galaxies belonging to different luminosity ranges or having different morphological types are 
different. These characteristics or 'marks' permit to understand the galaxy catalogs that carry all this information as realizations of 
marked point processes. Many attempts have been presented to quantify the dependence of the clustering of galaxies on their inner 
properties. 

Aims. The present paper summarizes methods on spatial marked statistics used in cosmology to disentangle luminosity, colour or 
morphological segregation and introduces a new one in this context, the mark connection function. 

Methods. The methods used here are the partial correlation functions, including the cross-correlation function, the normalised mark 
correlation function, the mark variogram and the mark connection function. All these methods are applied to a volume-limited sample 
drawn from the 2dFGRS, using the spectral type 77 as the mark. 

Results. We show the virtues of each method to provide information about the clustering properties of each population, the dependence 
of the clustering on the marks, the similarity of the marks as a function of the pair distances, and the way to characterise the spatial 
correlation between the marks. We demonstrate by means of these statistics that passive galaxies exhibit stronger spatial correlation 
than active galaxies at small scales (r < 20 h' 1 Mpc), and that the price galaxies have to pay to be close together is having smaller 
values of the assigned marks, which means in our case being more passive. Through the mark connection function we quantify the 
relative positioning of different types of galaxies within the overall clustering pattern. 

Conclusions. The different marked statistics provide different information about the clustering properties of each population. Different 
aspects of the segregation are encapsulated by each measure, being the new one introduced here -the mark connection function- 
particularly useful for understanding the spatial correlation between the marks. 

Key words. Cosmology: large-scale structure of Universe — Methods: data analysis — Methods: statistical 
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1. Introduction 

Galaxies of different morphological types show different clus- 
tering properties. It is well known, for example, that elliptical 
galaxies are preferentially found in high density environ ments, 
such as the centres of rich galaxy clusters (Dressier 1980), while 
the dominant population of the field are mainly spiral galaxies 
dDavis & GelleJl97dlDre"ssleill9 80). Second order characteris- 
tics as the two point correlation function have been used to quan- 
tify the clustering of galaxies with different morphologies, dif- 
ferent spectral characteristic s, different colours or belonging to 
differ ent luminosity ranges dPhillipps & Shanksll 19871; Ham ilton 
1988 ':lDavis et alJl988tLovedav et al.lll995LlHermit et alJ l996: 
iGuzzo et al.lll9971) . Bright galaxies show stronger spatial cor- 
relation than faint ones. Other clustering measures have been 
also used to quantify t he luminosity or morphological seg- 
regation: multifract als (Domfngu ez-Tenreiro & Martfnez| [l989; 
Domfngu ez-Tenreiro et al.l 1 19941). void probability functions 
(IVogelev et al.ll99lUCroton et al.l2004l). distributio ns of the dis- 
tances to the nearest neig hbours (lSalzer_et_alj[l990), etc. 

The two-point correlation function £(r) measures the excess 
probability of finding a neighbour at a distance r from a given 
galaxy when compared with that probability for a homogeneous 
Poisson process. Morphological segregation is encapsulated by 
the behaviour of £(r) when it is calculated separately for different 
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populations of galaxies. Elliptical galaxies show at small scales 
a correlation function with steepe r slopes and larger amplitudes 
than spirals (lLovedav et al.lfl995l) . A recent analysis of the Two 
Degree Field Galaxy Redshift Survey (2dFGRS) has shown the 
same trend when comparing populations for different spectral 
types, being the two-point correlation function steeper for pas - 
sive galaxies than for ac tive galaxies (Madgwi ck et al.l 12003). 
Also. IZehavi et alJ (120021) have analysed the distribution of red 
and blue galaxies in the Sloan Digital Sky Survey (SDSS) by 
means of the projected correlation funcion w p (r p ) showing that 
red galaxies display a more prominent and steeper real-space 
correlation function than blue galaxies do. 

The galaxy distribution can be considered a realisation of a 
point process. However, in many situations, each galaxy (point 
in the process) carries additional information regarding a given 
characteristic (e.g. morphological type) or a given numerical 
value that measures a given galaxy property: luminosity, colour, 
spectral type. If we attach this characteristic (mark) to the point 
in the process, we end up at a marke d point process, as it is 
called in mainstream s patial statis tics ( Stovan & Stovan|["l994l: 
iMartmez & Saarll2002l: llllian et al.ll2008~ 

In this work, we compare different statistical methods for 
the study of the marked galaxy distribution. We also introduce 
- for the first time in this context - the mark connection func- 
tion. We illustrate the usefulness of these methods by applying 
them to a volume-limited sample drawn from the 2dFGRS with 
marks given by the galaxy spectral type. In Section 2, we de- 



2 



V.J. Martinez et al.: Mark connection function 



scribe the sample and the marks assigned to the galaxies. In 
Section 3 we describe the different statistical methods consid- 
ered, and in Section 4 we show the results of applying them to 
our galaxy sample. In the conclusions, we stress the capabili- 
ties of the mark connection function to characterise the spatial 
correlation between the marks. 



2. The 2dFGRS subsamples 

To illustrate the different mark clustering measures, we used a 
nearly volume-limited sample drawn from the 2dFGRS and pre- 
pared by the 2dF team dCroto n et al. 2004). It contains galaxies 
with absolute magnitudes in the range -20 < Mb, < -19 at 
redshifts z < 0.13. In order to avoid the effects of complicated 
boundaries while using a simple estimator, we selected galaxies 
inside a rectangular parallelepiped inscribed in the North slice of 
2dFGRS. The final sample used contains N = 7741 galaxies and 
covers a volume of V ~ 10 6 (hr x Mpc) 3 where h is the Hubble 
constant in units of 100 km s _1 Mpc -1 . 

We characterized the galaxies in the sample using the spec- 
tral classification parameter rj (Madgwick et alj [20021) . Lower 
values of rj correspond to more passive or 'early-type' galax- 
ies, while larger values correspond to active or 'late-type' ones. 
In order to avoid negative values of the marks, we defined the 
mark used as m = rj + 10. This shift does not affect our con- 
clusions. Based on this 77 p arameter, we div i ded o ur sample in 
two populations, following Madg wick et al.l d2003l) : population 
'1' (passive galaxies) with 77 < -1.4, and population '2' (active 
galaxies) with rj > -1.4. These subsamples contain Ni = 3828 
and N2 = 3913 galaxies, respectively. We show the sample used 
in Fig.Q] 

In order to test the existence of mark segregation, we com- 
pared the results obtained for the different statistics with ran- 
dom relabelling simulations. In these, we keep the original po- 
sitions of galaxies, but redistribute the marks randomly among 
them. This corresponds to a model in which clustering is inde- 
pendent of the mark, or spectral type, of the galaxies. We simu- 
lated n = 200 realizations using the random relabelling method, 
and obtained their maximum and minimum values as function 
of distance r for each statistic. Deviations of the observed statis- 
tics from this range of values correspond to rejection of the 
mark-independent clustering model at a pointwise significance 
of 1 - ^ 99% (llllian et al.ll20"08h . 



3. Clustering analysis methods 

Recently, the clustering dependence on luminosity, colour or 
morphology has been analysed by means of the marked cluster- 
ing statistics, that allow to study the galaxy clustering as a func- 
tion of their properties, and moreover this approach provides us 
with different measures of the correlation betwee n the galaxy 
properties and the environment (Skibba et al. 2009). The galaxy 
distribution is interpreted as a realisation of a marked point pro- 
cess X M = {(Xi,rrii)}, where the mark m, denotes an intrinsic 
property of the galaxy at position x,. The mark can be the lumi- 
nosity, the spectral type, the colour, etc. In general, today galaxy 
catalogues provide quantitative marks ranging in a continuous 
interval rather than just a discrete characteristic like a galaxy be- 
ing spiral or elliptical. In any case, we shall also show how to use 
interesting second-order measures to disentangle clustering de- 
pendent characteristics of two populations by dividing the sam- 
ple in two parts using a significant value of the mark as threshold 
m t hres and separating the two populations according to the value 



of the mark: population '1' with m, < mthres an d population '2' 

With Mi > fWthres- 

We describe below the different methods we used to obtain 
information about galaxy clustering segregation. They are the 
classical partial correlation functions (for two discrete popula- 
tions), the normalized mark correlation function and the mark 
variogram (based on the use of continuous marks), and finally 
the mark connection function (based on the use of discrete 
marks). 

We computed the different statistics based on the estima- 
tion of the second-order intensity f unction for the unmarke d 
point process^ (Aajr)) pre sent ed in IStovan & Stovanl i 19941) . 
Pons-Borderfa et al] (Il999l) . and llllian et al.l (120081) . 
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where are the positions of the points, k(-) is a kernel function, 
and V(W n W r ) is the volume of the window (the parallelepiped 
in our case) intersect ed with a version of itself shifted by the 
vector r (see Fig. 1 in IPons-Bo rderfa et aTl l 19991) . 

In all our calculations we used the Epanechnikov kernel, 

fl - 4) for - w < x < 
otherwise 
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with width w = 1 h Mpc, and sampled the different functions 
with a step in r of 0.5 hr x Mpc. T his compact kernel is ver y 
well suited for correlation analysis (IPons-Borderfa et a"D [l999). 
We note, however, that the choice of a given kernel is not crucial, 
while the choice of the bandwith, w, is more important and plays 
the role of the binning in the standard calculation of correlation 
functions, where a top-hat kernel is typically used as default. 

3.1. Partial two-point correlation functions 

In the standard clustering analysis of the galaxy distribution, the 
two-point correlation function, £(r), measures the clustering in 
excess (£(r) > 0) or in defect (£(r) < 0) relative to a Poisson dis- 
tribution, for which ^(r) = 0. Whenever we want to compare the 
clustering properties of different populations of galaxies encap- 
sulated by their spatial correlations, we can consider the corre- 
lation function restricted to a given population, which are called 
partial correlation functions. In fact, for two populations of in- 
terest, one can consider three partial two-point correlation func- 
tions, namely £n(r), %2i(r), and ^n{r). The first two are those 
mentione d above for typ es 1 and 2, while the cross-correlation 
function dPeebleslll9 80) %vi(r) measures the excess probability 
of finding a neighbour of type '1' at distance r from a given 
galaxy of type '2', or vice versa. 

Based on equation (Q]), the partial two-point correlation func- 
tions were estimated as 
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1 Note that the relation between A 2 (r) and the standard correlation 
function used in cosmology is £(r) is A2(r) = n 2 [l +£(r)], where is the 
number density. The function g(r) = 1 + £(r) is known as the pair cor- 
relation function in spatial statistics. We use the convention of denoting 
the estimators by putting a hat'on top of the symbol of a given function 
to distinguish the estimator /ti(r) from the theoretically defined func- 
tion A 2 (r). Although this is not standard in cosmology, it is an extended 
convention in spatial statistics, and it is quite usefu l when different esti- 
mator s of a single function are discussed (see, e.g.,|Pons-Borderi a et al.1 
C99l)). 



Fig. 1. Tridimensional plot of the galaxy sample used. Red dots correspond to early-type galaxies (population '1'), and blue dots to 
late-type galaxies (population '2'). The parallelepiped dimensions are 254 x 133 x 31 hr x Mpc. 



where xP are the positions of galaxies of population i, and «, = 
NJV. 

We estimated the error of the measured correla tion functions 
using the jackknife method (Norberg et al. 2009). We divided 
the data volume in 32 equal, nearly cubic, sub-volumes. We gen- 
erated the corresponding 'mock' datasets omitting one of these 
sub-volumes at a time, and calculated the correlation functions 
for these. The jackknife errors for each scale, cr (/ (r), are then 
obtained as 



2, v N k -l 



Nt 



2(4(r)-| y (r)) , 
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where £fXr) is the partial correlation function ^ij(r) of the 'mock' 
dataset k, £,y(r) is the value averaged over these datasets, and 
N k = 32. 



3.2. Normalized mark correlation function 

Stov an & Stovanl d 1994b introduced the normalized mark cor- 
relation function. To define this function let us first define the 
quantity 



A^Kx], mi), (X2, m^)\dV\dm\dVidm,2 



(3) 



as the joint probability that in the volume element dV\ lies 
a galaxy with the mark in the range [m\,m\ + dm{\ and that 
another galaxy lies in dV2 with the mark in [ni2, ni2 + dni2] 
(Mart inez & Sa ar 2002). The normalized mark correlation func- 



tion is 



t 2 A 2 (r) 
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mim2/i^((jti, mi), (X2, mii)dm\dmi, (4) 



for A%(r) + 0, where in is the mean of the marks. 

Despite its name the mark correl ation function is not a strict 
correlation function (Schlather 2001), but it describes important 



aspects of the spatial correlations of marks. A true mark corre- 
lation is a function given by Eq. but replacing the product 
m\m.2 by the product of the differences (ni\ - m)( m 2 - m). The 
normalizing denominator m 2 must then be replaced by <x 2 , the 
variance of the marks. In any case, k mm (r) < 1 represents inhi- 
bition of the marks at the scale r. For example, in forests it is 
typically found that trees with larger stem diameter (mark) tend 
to be isolated, since they make use of much more ground and 
sun-light resources than smaller trees. Using luminosity as the 
mark, the opposite effect has been f ound for the galaxy distribu - 
tion, i.e., k mm (r) > 1 at small scales (Beisba rt & Kerscherl2 000). 
implying stronger clustering of brighter galaxies at small sepa- 
rations, in a greement with p revious results showing this kind of 
segregation (Hamilton 1981). 

We estimated the normalized mark correlation function as 
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3.3. Mark variogram 

The mark variogram, y{r) (IWaider & Stovanl [l996l: 
Beisb art & Kerscheri l2000l) . is a measure of the similarity 
of the marks depending on the distance between galaxies. It is 
defined as 
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2A 2 (r) 



II (m < 



-m2) 2 A^(( Xl . 



mi), (X2, ni2))dm\dni2 .(6) 



When the clustering properties of a marked point pattern are in- 
dependent of the marks, then the mark variogram y(r) is con- 
stant and takes, naturally, the value of the variance, (T 2 m , of the 
mark distribution. In the presence of segregation, the fact that 
y(r) > tj 2 n indicates that galaxy pairs at distance r tend to have 
different marks, while the contrary, y(r) < cr^, is an indication 
that these galaxy pairs tend to have similar marks. 
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We estimated the mark variogram as 



f(r) 



1 ^ V ( m i- m j) 2k ( r ~ l*> 



V(W n W Xi - Xi ) 



(7) 



3.4. Mark connection function 

A statistical tool to characterize the spatial correlation between 
the marks of a point pattern with discrete marks is the mark con- 
nection function pij(r), which represents the conditional prob- 
ability to find two galaxies of type i and j at positions sepa- 
rated by a distance r, under the condition that at these positions 
there are indeed galaxies. This function yields information dif- 
ferent to that from the partial correlation functions, £,y(r), as 
shown, for example, in llllian et akl (120081) . By its definition it 
gives the relative frequencies of mark pairs (i, f) of distance r. 
While Hijif) takes large values if there are many (/, y')-pairs at 
distance r, Pij(r) is large if the proportion of (i, y)-pairs in all 
pairs at distance r is large. So it may happen that for some r, 
£ij(r) has a minimum, but Pij(r) has a maximum, if there is only 
a small number of point pairs at distance r in the whole pattern, 
but many of them are exactly (i, j')-pairs. Experience shows that 
often Pij(r) is able to find finer structures in point patterns than 
^ij(r), because of the nature of Pij(r) as a conditional probability. 

If the marking is independent of clustering, then Pij(r) are 
constant, 



(8) 



Here p-, is the probability that a randomly chosen galaxy is of 
type i. The p, are estimated as 

We calculated Pij(r) based on the estimation of the partial corre- 
lation functions as 



Pij(r) = pipj 



W + 1 

|(r) + 1 



where £(r) is the two-point correlation function of the full sam- 
ple. 



4. Results 

4.1. Partial two-point correlation functions 

Fig. [2] shows the three corresponding partial two-point correla- 
tion functions, estimated according to Eq. ©. All three show 
clearly the high degree of clustering within the pattern of galax- 
ies. It is obvious that the correlation function for the type '1' 
passive galaxies is steeper than for the type '2' active galaxies as 
well as for the (1,2) pai rs. This result c orroborates the spectral 
segregation detected by Madgwic k et al.l d2003l) for the 2dFGRS. 

4.2. The normalized mark correlation function 

The k mm (r) for our sample, estimated according to Eq. (|5j, is 
shown in Fig. [3] The curve for k mm (r) shows a weak negative 
correlation or spatial inhibition: k mm {r) < 1. The range of corre- 
lation is about 20 h~ l Mpc, where k mm (r) gets values close to 1. 
It is interesting to compare this result with the £ mm (r)-function 
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Fig. 2. The partial two-point correlation functions, ^y(r), esti- 
mated for population '1' (early-type) and population '2' (late- 
type) galaxies in our sample. The top panel shows the three func- 
tions together. The three lower panels show each of them sepa- 
rately (solid lines with error bars estimated using the jackknife 
method), together with a shaded band showing the minimum and 
maximum values for the 200 realizations of the random rela- 
belling simulation. The dot-dashed lines correspond to £(r) for 
the full sample, which is the expected value of all £y(r) in the 
absence of segregation. 



shown in Beis bart & Kerscherl (l2000b using the galaxy absolute 
luminosity L as the mark. They obtain an increasing behaviour of 
kmm(r) at small scales with k mm (r) > 1 for r < 12h~ l Mpc, show- 
ing that bright galaxies are stronger correlated than faint ones. In 
our case, the tendency of the values of k mm {r) to be smaller than 
1 at short scales indicates that the price galaxies have to pay for 
being close together is to have reduced values of the marks, i.e., 
being more passive. 

4.3. The mark variogram 

In Fig.H] we show the mark variogram for our sample, obtained 
according to Eq. (|7J. This function is monotonously increasing. 
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Fig. 3. Normalized mark correlation function, k m „,(r), for our 
sample (solid line). The shaded band shows the minimum and 
maximum values for the 200 realizations of the random rela- 
belling simulation, while the dot-dashed line correspond to the 
value for the case with no segregation, k mm (r) = 1. 

In this case the interpretation is straightforward: y(r) shows that, 
for separations r < 10 h~ x Mpc, galaxy pairs tend to have similar 
marks, that is, similar spectral type. 

This result is partially explained with the previous one shown 
by the k mm function: galaxies close together exhibit smaller val- 
ues of the attached mark (spectral type). 
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Fig. 4. Mark variogram, y(r), for our sample (solid line). The 
shaded band shows the minimum and maximum values for the 
200 realizations of the random relabelling simulation, while the 
dot-dashed line corresponds to the value for the case with no 
segregation, y(r) = a 3 ^ - 6.25. 



4.4. The mark connection function 

We show the Pij(r) obtained for the 2dFGRS galaxies, together 
with the results of our random relabelling simulations, in Fig. 
The first panel shows very neatly that, for scales r < 20 h~ l Mpc, 
the clustering of early-type galaxies is stronger than the clus- 
tering of late-type galaxies. The three bottom panels show that 



the deviation of the observed py(r) from the case of random la- 
belling is significant at these scales. 

Moreover, the figure shows clear differences in the spatial 
correlations of galaxies of the two types. In an overall clustering 
of all galaxies, we can outline that: 

1 . Galaxies of type ' 1 ' (passive or early-type) are strongly clus- 
tered up to distances of 20 h~ x Mpc. 

2. The conditional probability to find two galaxies of type '2' 
(active or late-type) at two positions separated by a distance 
r (under the condition that at these locations are galaxies) is 
smaller than the same probability for random labelling of the 
marks for scales r < 20 hr x Mpc. 

3. Galaxy pairs having one member of type '1' and the other 
member of type '2' are less frequent than for random la- 
belling up to distances of 10 hr x Mpc. 

In summary, all galaxies form a highly clustered pattern. In 
this pattern, the passive galaxies tend to be close to other passive 
galaxies, while positioning of active galaxies is less affected by 
other active galaxies. However, they tend to avoid positions close 
to passive galaxies. 

This shows clearly the power of the mark connection func- 
tion as an analytical tool in comparison to the partial pair corre- 
lation function. While, for the untrained eye, the curves in Fig. 
2 are quite similar and show little structure, the curves in Fig. 5 
give valuable information about the inner structure of the mark 
distribution. Obviously, the idea to consider characteristics of 
the nature of conditional probabilities helps to divulge structural 
details which would be otherwise overlooked. 

The problem are the mutual positions, given the positions 
of all galaxies without mark information. Since the three partial 
two-point correlation functions, shown in Fig. 2, are different 
for a large range of scales, the marking with marks 1 and 2 can 
not be an independent marking, where every galaxy obtains its 
mark randomly, independent of the other galaxies. In contrast, 
there must exist a spatial correlation between the marks. As it 
has been shown in Fig. 5, the mark connection function is the 
appropriate tool to measure this correlation. 

5. Conclusions 

We have used a volume-limited galaxy sample from the 2dF- 
GRS to test different statistical measures used to disentangle 
mark segregation in the distribution of the galaxies. The mark 
attached to each galaxy of the sample was its spectral type r\. 
For some of the statistics, the value of the mark enters directly 
into the functions used for measuring segregation: the normal- 
ized mark correlation function k mm (r) and the mark variogram 
y(f). For other functions, as the partial correlation functions or 
the mark connection function, the sample has been split into two 
populations corresponding to passive or early-type galaxies with 
i] < -1.4 and active or late-type galaxies with rj > -1.4. We 
summarise our results as follows: 

1. The partial correlation functions, including the cross- 
correlation function, inform us about the degree of clustering 
of each population separately. It shows that passive galaxies 
exhibit stronger clustering at small separation. Nevertheless, 
there is no information about the spatial correlation between 
the marks. 

2. The normalized mark correlation function shows that having 
smaller values of the marks, i.e., smaller values of spectral 
type (being more passive), is a clear condition for galaxies 
being close to each other in the overall clustering pattern. 
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Fig. 5. Mark connection functions, pij(r) obtained for 'early- 
type' (population '1') and 'late-type' (population '2') galaxies 
in our sample. The top panel shows the three functions together. 
The three bottom panels show pu(r), P22W, and pnir) sepa- 
rately (solid lines), together with the shaded band showing the 
minimum and maximum values for the 200 realizations of the 
random relabelling simulation. The dot-dashed lines correspond 
to the expected values for the random labels case, according to 
Eq.[8l 



Applied on the 2dFGRS volume-limited sample, the mark 
connection function clearly shows that passive galaxies are 
clustered up to distances of 20 h^ 1 Mpc, while active galax- 
ies exhibit weak spatial anticorrelation of the mark up to dis- 
tances of 20 h Mpc. Mixed pairs are less frequent up to 
distances of 10 h' 1 Mpc. 
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3. The mark variogram, in addition, shows that at small separa- 
tions galaxy pairs tend to have similar marks. 

4. The mark connection function has been introduced here for 
the first time in the analysis of the marked galaxy distribu- 
tion. The function measures the conditional probability 
to find at two positions, separated by a distance r, a galaxy 
of type 'i' and a galaxy of type 'j' under the condition that 
at these positions there are indeed galaxies. This function 
yields information different from that of the partial correla- 
tion functions £y(r). This more sophisticated measure, hav- 
ing a nature of conditional quantities, is an efficient statis- 
tical tool to characterize the spatial correlation between the 
marks, filtering out the relative frequencies of the mark pairs 
(i, j) at distance r. 



